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WIND-TUNNEL INVESTIGATION OF THE EFFECTS OF WING THICKNESS 
ON THE STATIC LONGITUDINAL AND LATERAL STABILITY 
OF UNSWEPT WINGS OF ASPECT RATIO 3 
AT HIGH SUBSONIC SPEEDS 

By William C. Hayes, Jr., and Edward C. Polhamus 


SUMMARY 


An investigation was made to determine the effects of wing-thickness 
on the static longitudinal and lateral stability of a wing-fuselage com- 
bination at high subsonic speeds. The wings had an aspect ratio of 3, 
a taper ratio of 0.5> zero sweep of the half-chord line, and NACA 65A002, 
NACA 65A004, and NACA 65AOO6 airfoil sections parallel to the plane of 
symmetry. The variations of side force, yawing moment, and rolling 
moment with angle of sideslip were obtained as a continuous trace through 
the angle-of-sideslip range in order to detect any nonlineaxities. 

Below a Mach number of 0.90, wing thickness had little effect on the 
lift-curve slope; however, above a Mach number of 0.90 a considerable 
increase was obtained for the 2-percent-thick wing. An increase in max- 
imum lift with decrease in wing thickness was noted throughout the Mach 
number range. A rearward movement of the aerodynamic -center location 
occurred at a lower Mach number than that predicted by theory. A decrease 
in drag due to lift with increase in wing' thickness was noted through the 
Mach number range to O.92, whereas at a Mach number of 0.95 the reverse 
Was true . 

The values of the rate of change of the effective dihedral parameter 
with lift coefficient were in fair agreement with theory except for a 
moderate increase with Mach number which was not predicted and which was 
greatest for the 6-percent-thick wing. Although wing thickness had little 
effect on directional stability, all three wing-fuselage combinations 
provided a desirable reduction in the wing-fuselage directional insta- 
bility at the higher angles of attack. 
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INTRODUCTION 


An extensive research program is being conducted in the Langley 
high-speed 7- by 10-foot tunnel to determine the effects of wing geometry 
on the static stability characteristics of wing-body combinations at 
high subsonic speeds. This program has included a systematic investiga- 
tion of the effect of sweepback (ref. l), aspect ratio (ref. 2), taper 
ratio (ref. 3), and geometric dihedral (ref. 4) on the longitudinal and 
lateral stability characteristics for Mach numbers up to about 0.95. 

In addition, the static stability of two delta plan-form wings has been 
investigated (refs. 5 and 6). A summary of the rolling moment due to 
sideslip characteristics for all these investigations has been published 
in reference J. 

Although the program is fairly complete with respect to plan-form 
characteristics, the effects of wing thickness have not been studied. 

In addition, the angle— of— attack range was rather limited — especially 
at the higher Mach numbers - and since (except for some limited data at 
low angles of attack) the lateral data were obtained at fixed sideslip 
angles (p = ±4°), no indication of possible nonlinearities in directional 
stability in the very low sideslip range was obtained. Therefore, the 
purpose of the present investigation was to determine the effect of wing 
thickness on the static stability characteristics through a large angle - 
of -attack range for Mach numbers up to 0.95. Testing through a large 
angle -of -attack range was made possible by use of considerably smaller 
models than those used in the previous investigations; thereby, the 
tunnel power and choking limitations were alleviated. Also, in order 
to study possible nonlinearities in the low sideslip r ang e, the models 
were tested on a variable-yaw sting which was designed to allow contin- 
uous records of the forces and moments to be obtained through the side- 
slip range. ° 

Three wings having an aspect ratio of 3, a taper ratio of 0.5, zero 
sweep of the half- chord line, and NACA 65A002, NACA 65A004, and NACA 65AOO6 
airfoil sections parallel to the plane of symmetry were tested. 


COEFFICIENTS AND SYMBOLS 


In figure 1 are shown the axis system and the positive directions 
of forces, moments, and angles used in presenting the data. With the 
exception of lift and drag, the data are presented about the body axes 
of the model. The center of moments is located at 25 percent of the 
mean aerodynamic chord of the wing. 
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Cl 


Cd 


ACd 


Cm 

Cz 

Cn 

C Y 

<1 

P 

Vo 

M 

S 

b 


c 


lift coefficient, 

qS 

drag coefficient, 

qS 

total drag coefficient minus drag coefficient at zero lift 

pitching-moment coefficient, g^cM.ng ™^nt 

qSc 

rolling-moment coefficient, moment 

qSb 

yawing-moment coefficient, Yawing moment 

qSb 

lateral-force coefficient, L ateral force 

qS 

PV ^ 

dynamic pressure, — lb/sq ft 

mass density of air, slugs/cu ft 
free-stream velocity, ft/sec 

Mach number 


wing area, sq ft 
wing span, ft 


wing mean aerodynamic chord. 



c wing chord at any spanwise station, ft 

t maximum wing thickness at any spanwise station, ft 

a angle of attack, deg 

0 angle of sideslip, deg 

A c /2 angle of sweepback of half -chord line, deg 

R Reynolds number 
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A 

A 


wing aspect ratio, b^/s 
wing taper ratio, 

Root chord 


X ac distance to aerodynamic center measured from leading edge of 

wing mean aerodynamic chord 

R . 

ba, 

C Yp 


Cn. 


3 


' l P 


lateral-force coefficient due to sideslip, per deg 

dp 

yawing-moment coefficient due to sideslip, per deg 

dp 

rolling-moment coefficient due to sideslip, per 


rate of change of with lift coefficient, dCL 


MODEL AND APPARATUS 


Details of the complete model are presented in figure 2. The fuse- 
lage consisted of a cylindrical aluminum center section fitted with 
plastic nose and tail sections. With this fuselage were tested three 
stainless-rsteel wings having an aspect ratio of 3 , a taper ratio of 0 . 5 , 
zero sweep at the half chord, and NACA 65A002, NACA 65A004, and 
NACA 65 AOO 6 airfoil sections parallel to the plane of symmetry. 

Schematic diagrams and photographs which illustrate the mechanical 
system used to obtain the continuous yaw records are presented in fig- 
ure 3 . With this system an angle -of -side slip range of approximately 
- 5 ° to 5 ° and an angle -of -attack range of - 3 ° to 23 ° can be obtained by 
remote operation. The prime motivation of the system is a l/4-horsepower 
electric motor which drives the piston within the master cylinder. The 
amount and direction of the movement of this piston is automatically 
controlled by a system of limit switches and clutches and is duplicated 
by the piston within the slave cylinder. Longitudinal movement trans- 
mitted from the slave cylinder to the sliding block is translated to 
lateral motion by means of a diagonal slot in the sliding block and a 
pivoted follower arm. The rate of change in sideslip angle p is 
approximately 1 degree per second which is believed to be slow enough 
to eliminate any appreciable effects of transient phenomena. Forces and 
moments are measured with a six- component strain-gage balance mounted 
within the model. 
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TESTS AND CORRECTIONS 


Tests were made in the Langley high-speed 7- by 10-foot tunnel for 
a Mach number range from 0.60 to 0.95 corresponding to a Reynolds number 

range of 1.15 x 10 ^ to 1.50 x 10 ^ based on the wing mean-aerodynamic- 
chord length (0.299 foot). An angle -of -attack range of -5° to 23° was 
obtained at the lower Mach numbers. At the higher Mach numbers the 
angle-of -attack range was limited (0° to approximately 13° at M = O. 95 ) 
to avoid exceeding the design- load limits of the balance. An angle-of- 
sideslip range of approximately - 5 ° to 5 ° was obtained at all test angles 
of attack. Jet -boundary corrections determined by the method of refer- 
ence 8 and blocking corrections determined by the method of reference 9 
were found to be negligible and therefore were not applied to the data. 
The angle of attack and angle of sideslip have been corrected for deflec- 
tion of the sting support system and balance under load. 


RESULTS AND DISCUSSION 
Longitudinal Stability Characteristics 


The basic lift, drag, and pitching -moment characteristics of the 
fuselage alone and the wing-fuselage combinations are presented in fig- 
ures 4 to 8 . The unrealistic drag coefficients indicated for the fuse- 
lage alone at low angles of attack are due to the fact that the drag on 
the fuselage alone was of the order of the accuracy of the balance which 
was designed to measure the large forces encountered by the wing-fuselage 
combination at high angles of attack. The upper part of figure 9 pre- 
sents, as a function of Mach number, a comparison of the experimental 
lift-curve slopes for the three different wing thicknesses with esti- 
mates based on thin-airfoil theory. The variation of the theoretical 
lift-curve slope with Mach number was calculated by the method presented 
in appendix A of reference 7 using the method of reference 10 to account 
for the wing-fuselage interference. It will be noted that below a Mach 
number of about 0.90 there is relatively little effect of wing thickness 
and that the theory is in fair agreement with experimental values; how- 
ever, above a Mach number of 0.90 the 2-percent-thick wing exhibited a 
considerably more rapid increase in lift-curve slope with Mach number 
than do the ^-percent- and 6 -percent-thick wings. This is apparently 
due to the ability of the thinner wing to maintain (at an angle of 
attack) supersonic flow over a larger portion of the chord. (See 
refs. 11 and 12 .) 

The maximum lift values presented in the lower part of figure 9 
represent the first peak and neglect the recovery of slope at the higher 
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angles of attack. The results indicate a rather large increase in max- 
imum lift with decreasing wing thickness throughout the Mach number range 
investigated. It will be noted, however, that the variation with Mach 
number for the 2-percent-thick wing is somewhat different from that for 
the 4-percent- and 6-percent-thick wings, with a decrease occurring 
between 0.6 and 0.8. All three wings exhibit a rather rapid increase 
in maximum lift at the higher Mach n umb ers - 

Figure 10 presents the variation of the aerodynamic-center location 
(aft of leading edge pf mean aerodynamic chord) with Mach number for the 
three wings investigated. Also shown are the theoretical thin-airfoil 
estimates which were determined by the method of reference 10. (The method 
of reference 10 was modified slightly with regard to the subsonic wing 
alone center of pressure presented in figure 17 of reference 10 since the 
values differ considerably from those determined by the lifting- surface 
theory. Values from references 13 and l4 were used to replace the 
A = 1 curve and were used as a guide to fair the curves for the other 
taper ratios into the slender -body values at |3A = 0.) The experimental 
results indicate general agreement with theory in that the forward move- 
ment of the aerodynamic center associated with increasing induced camber 
is followed at the high Mach numbers by a rapid rearward movement asso- 
ciated with the establishment of chordwise loading of the supersonic type. 
It will be noted, however, that the experimental results for all three 
thickness ratios indicate a considerably more rapid forward movement and 
an earlier rearward movement than indicated by the theory. Both of these 
deviations are probably associated with the fact that the thickness- and 
lift-induced perturbation velocities cause an early establishment of local 
supersonic velocities on the wings which is neglected in the theory. The 
lack of progression toward the theory as the wing thickness is reduced is 
probably due to the fact that a reduction in thickness is accompanied by 
an increase in the effectiveness of lift coefficient in reducing the criti- 
cal Mach number. This effect can, even at the low lift-coefficient range 
over which the aerodynamic-center locations were determined, cancel the 
favorable thickness effect which exists at zero lift. (See ref. 15 . ) 

The more rearward location of the aerodynamic center of the 2-percent- and 
4-percent-thick wings relative to that for the 6-percent-thick wing may be 
associated with leading-edge separation which occurs on thin airfoils. 

Figure 11 presents the variation of the drag-due -to- lift parameter 
with Mach number for each of the experimental wings compared in each case 
with the theoretical values for full leading-edge suction l/itA and zero 
leading-edge suction l/C]^ (see ref. l6), where Ci^ was determined 

from data of the present investigation. It will be noted that at Mach 
numbers up to approximately 0-92 the drag due to lift increased with 
decrease in wing thickness probably because of loss of leading-edge suc- 
tion associated with the decrease in leading-edge radius. At the highest 
test Mach number the drag due to lift decreased with decrease in wing 
thickness because at this Mach number the resultant force i-s normal to 
the wing chord and, therefore, the thin wing with a higher lift-curve 
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slope has less drag due to lift. It is noted that in the case of the 
2 -percent-thick wing, the experimental values of drag due to lift are 
higher than the theoretical values which assume zero leading-edge suction. 
This is so probably because of the increase in profile drag associated 
with the forward movement of the transition with increase in lift 
coefficient. (See ref. 17 .) 


Lateral Stability Characteristics 

The basic lateral-force, yawing -moment, and rolling-moment charac- 
teristics of the fuselage alone and the wing-fuselage combinations are 
presented in figures 12 to 15 . As mentioned previously, the lateral- 
stability data were determined from continuous records of the force and 
moment variation with sideslip angle. The static-lateral-stability 
derivatives C Y p, C n p, and C^ were computed at each angle of attack 

by taking the algebraic difference between the values of the particular 
component at each of the limiting angles of sideslip (0 « ± 5 °). These 
values were then divided by the algebraic difference between the angles 
of sideslip. This is illustrated in sketch (a) for a linear variation of 



Sketch (a) 
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In some cases, however, the variations of side force, yawing moment, 
and rolling moment were nonlinear with angle of sideslip. Sketch (b) 
shows two sample traces of nonlinear variations of rolling moment with 
angle of sideslip. It should be noted that the same method was used to 
determine the derivatives as was employed in the cases of linear varia- 
tion. 



Sketch (b) 


In order to differentiate between the derivatives obtained from 
linear or nonlinear variations with p, different types of symbols are 
used in figures 13 to 15* In these figures, open symbols are used to 
denote the angles of attack at which the variation of force or moment 
with angle of sideslip was considered linear (sketch (a)), whereas the 
solid symbols denote those angles of attack at which the variation was 
considered nonlinear (sketch (b)). 


Inasmuch as the effect of wing thickness appears to be greatest for 
the effective dihedral parameter C ^ and since the directional stability 

and side-force characteristics (c n p and Cyp> respectively) will be 
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affected to a large extent by the particular tail assembly used, the 
present discussion is mainly concerned with the ef fective ' dihedral param- 
eter'. In view of the fact that Cjp is more closely associated with 

lift coefficient than angle of attack,, a comparison of the variation of 
Cjp with Cl for the different wing thicknesses is presented in fig- 
ure 16 for the various Mach numbers investigated. With the exception of 
the highest Mach number the results indicate the effective dihedral param- 
eter Cjp to be a fairly linear function of lift coefficient below a 

lift coefficient of about 0.40. It can also be noted that in this range 
the effect of wing thickness and Mach number is relatively small. This 
is more clearly indicated in figure 17 where the slope 

curve of C^ plotted against Cl is presented as a function of Mach 

number for the three wing thicknesses. Also shown is an estimate by the 
method presented in reference 7* The experimental values are in fairly 
good agreement with the estimated values at low Mach numbers; however, 
the experimental results indicate a moderate increase with Mach number 
which is not predicted by the method of reference 7- Although there 


appears to be essentially no effect of wing thickness on C; 


between 


the 2 -percent -thick and 4-percent-thick wings, the values for the 
6 -percent- thick wing are slightly larger throughout the Mach number range. 
It can be noted that at the higher Mach numbers where supersonic flow is 
fairly well established on the wings there is a decrease in as 


indicated in reference 7 for supersonic speeds. 


Regarding the effective dihedral characteristics at high lift coef- 
ficients the results (fig. l 6 ) indicate rather large effects of wing 
thickness, Mach number, and lift coefficient. The curves at the higher 
lift coefficients are characterized by a rapid increase in C 2 ^ with 

increasing lift coefficient. Wing thickness had a pronounced effect, 
with the increase occurring at higher lift coefficients as the thickness 
decreased. As the Mach number increased beyond 0.90, the effect of lift 
coefficient decreased, and at a Mach number of 0-95 little effect of lift 
coefficient or wing thickness occurred. 


As mentioned previously, the variation of rolling-moment coefficient 
C 1 with sideslip angle p was considerably nonlinear in some cases 
which are represented by solid symbols in the basic data. (See fig. 15 . ) 
This is illustrated in more detail in figure 18 where the characteristics 
of the 6 -percent- and 2-percent -thick wings are compared at a Mach number 
of O. 85 . For orientation with regard to lift coefficient, the upper part 
of the lift curve is also presented. Also shown are the continuous traces 
of rolling-moment coefficient as a function of sideslip angle from which 
the values of C^p were determined. From these traces it can be noted 
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that in the angle-of -attack range where there are large reductions in 
Cjp associated with the thin wing, the variation with sideslip angle 

is very nonlinear and the derivative appears to be of little value. A 
more complete set of traces for all three wing thicknesses at this* same 
Mach number is presented in figure 19 . 


Figure 20 presents the combinations of lift coefficient and Mach 
number for which nonlinear variations of rolling-moment coefficient with 
sideslip angle exist. 

In connection with the directional-stability parameter C np , it can 

be noted (fig. l4) that in addition to the fact that nonlinear variations 
with sideslip exist, the wing-fuselage instability generally decreases 
considerably at the higher angles of attack. This trend, which is desir- 
able in view of the usual loss of vertical-tail contribution at high 
angles of attack, is opposite that obtained for sweptback wings and 
appears to be associated with a favorable wing interference on the fuse- 
lage afterbody. (See ref. 18 . ) 


CONCLUSIONS 


The results of a wind-tunnel investigation to determine the effects 
of wing thickness on static longitudinal and lateral stability charac- 
teristics of unswept wings of aspect ratio 3 at high subsonic speeds 
indicate the following conclusions: 

1. Below a Mach number of O .90 the lift- curve slope is little 
affected by wing thickness and shows fairly good agreement with theory 
Above a Mach number of 0.90 the lift-curve slope for the 2-percent-thick 
wing exhibits a more rapid rise with increase in Mach number than do 
the thicker wings. Throughout the Mach number range the maximum lift 
decreases with an increase in wing thickness. 


2. The aerodynamic centers of the 2-percent- and ^-percent-thick 
wings are somewhat rearward of that of the 6 -percent-thick wing at sub- 
critical speeds. All three wings exhibit the trends with Mach number 
indicated by theory but these trends occur at considerably lower Mach 
numbers than theory indicates. 


,, The <3rag due to lift increases with decrease in wing thickness 

through the Mach number range from O .60 to approximately O. 92 . However 
at a Mach number of approximately 0-95 the drag due to lift decreases ' 
with a decrease in wing thickness. 

. ^ lues of the rate of change of the effective dihedral param- 

eter with lift coefficient are in fair agreement with theory except for 
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a moderate increase with Mach number which was not predicted. The 
6-percent-thick wing has a value which is slightly higher than those for 
the thinner wings. 

5. Wing thickness has relatively little effect on the directional 
stability; all three wing-fuselage combinations provide a desirable reduc- 
tion in the wing-fuselage instability at the higher angles of attack. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., May l4, 1956. 
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Figure 2.- Details of complete model. All dimensions are in inches 
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Figure 3.- Concluded. 
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Figure 5.- Variation of angle of attack with lift coefficient for wing 
fuselage combinations. Solid symbols indicate origins of axes. 
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Lift coefficient, C, Lift coefficient ,C L 
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Lift coefficient ,C L Lift coefficient ,C L 
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Figure 6.- Variation of drag coefficient with lift coefficient for wing- 
fuselage combinations. Solid symbols indicate origins of axes. 
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(a) t/c = 0.02. 

Figure 7 .- Variation of pitching -moment coefficient with angle of attack 
for wing-fuselage combinations . Solid symbols indicate origins of axes 
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Figure 8.- Variation of pitching -moment coefficient with lift coefficient 
for wing-fuselage combinations. Solid symbols indicate origins of axe, 
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Figure 10.- Comparison of theory with experimental variation of 

aerodynamic 'center . 
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Figure 13.- Variation of Cy^ with a for wing-fuselage combinations. 
Solid symbols indicate nonlinear variations of Cy with 3 . 
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Figure l4 . - Variation of Cn^ with a for wing -fuselage combinations 
Solid symbols indicate nonlinear variations of C„ with B. 
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Figure 15.- Variation of C lp with a for wing-fuselage combinations. 
Solid symbols indicate nonlinear variations of Cl with 0. 



NACA RM L56E30a 


CONFIDENTIAL 


39 



CONFIDENTIAL 


Figure 16 .- Variation of effective dihedral parameter with lift coefficient. 
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Figure 18.- Effect of wing thickness and angle of attack 
on Cip at M = O. 85 . 
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Figure 20.- Summary figure showing test limits and areas in which the 
variation of C i with (B was generally nonlinear. 



CONFIDENTIAL 


CONFIDENTIAL 


